Multi-impurity Anderson model for quantum dots coupled in parallel 
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The system of several (N) quantum dots coupled in parallel to the same single-mode conduction 
channel can be modelled as a single-channel A^-impurity Anderson model. Using the generalized 
Schrieffer- Wolff transformation we show that near the particle-hole symmetric point, the effective 
Hamiltonian in the local moment regime is the A''- impurity S — 1/2 Kondo model. The conduction- 
band-mediated RKKY exchange interaction between the dots is ferromagnetic and at intermediate 
temperatures locks the moments into a maximal spin S = N/2 ground state. We provide an 
analytical estimate for the RKKY interaction. At low temperatures the spin is partially screened by 
the conduction electrons to N/2 — 1/2 due to the Kondo effect. By comparing accurate numerical 
renormalization group results for magnetic susceptibility of the A'^— impuriy Anderson model to the 
exact Bethe-Ansatz results of a S = A'^/2 SU{2) Kondo system we show, that at low-temperature 
the quantum dots can be described by the effective 5* = A'^/2 Kondo model. Moreover, the Kondo 
temperature is independent of the number of impurities A'^. We demonstrate the robustness of the 
spin A'^/2 ground state as well as of the associated S — N/2 Kondo effect by studying the stability of 
the system with respect to various experimentally relevant perturbations. We finally explore various 
quantum phase transitions driven by these perturbations. 

PACS numbers: 72.15.Qm, 73.23.Hk 73.63.Kv, 71.10.Hf 
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I. INTRODUCTION 

The Kondo effect emerges as the increased scattering 
rate of the conduction band electrons at low tempera- 
tures due to the presence of magnetic impurities which 
induce spin-flip scattering. It leads to various anoma- 
lies in the thermodynamic and transport properties of 
the Kondo systems. It is usually described using simpli- 
fied quantum impurity models such as the Kondo model 
and the Anderson model^.The quantum impurity mod- 
els attract the interest of the solid state physics com- 
munity both due to their unexpectedly complex behav- 
ior and intrinsic beauty, as well as due to their ubiq- 
uitous applicability to a vast array of physical systems 
such as bulk Kondo systems, heavy- fermion compounds 
and other strongly correlated systemsSi, dissipative two- 
level systemsi*^, single magnetic impurities and quantum 
dotslS^. 

After the properties of single-impurity models were un- 
raveled using a complementary set of techniques (the 
scaling approach, Wilson's numerical renormalisation 
group, Bethe-Ansatz solution and various large- A'^ expan- 
sion schemes)^, the attention has increasingly focused to 
multi- impurity models. Research in this field has recently 
increased due to a multitude of experimental results 
made possible by advances in micro- and nanotechnol- 
ogy. The multi-impurity magnetic nanostructures under 
study are predominantly of two kinds: clusters of mag- 
netic adsorbates on surfaces of noble metals (Ni dimers^, 
Ce trimers^, molecular complexes^) and systems of mul- 
tiple quantum dotsiSiiLi^iiMl. 

The most important additional element that emerges 
in multi-impurity models is the Ruderman-Kittel- 
Kasuya-Yosida (RKKY) exchange interaction-- . It arises 



when the magnetic moments on the impurities induce 
spin polarization in the conductance band which leads 
to magnetic coupling of moments that are separated in 
space. The RKKY interaction depends on the inter- 
impurity distance and can be either ferromagnetic or an- 
tiferromagnetic. 

The competition between the antiferromagnetic 
RKKY interaction and the Kondo effect in two magnet- 
ically coupled local moments leads to a quantum phase 
transition at J '^ Tk between strongly bound local mag- 
netic singlet for J ^ Tk and two separate Kondo sin- 
glets for J < J^^i6,i7,i8,i9,20_ rpj^g ^^jg ^j ^^g antifer- 
romagnetic exchange interaction was also studied in the 
context of double quantum dots (DQD)Sii22i2iS12^. Two 
mechanisms can contribute to the effective exchange in- 
teraction between the dots: the conduction-band medi- 
ated RKKY interaction and the super-exchange mecha- 
nism due to inter-dot electron hopping. Depending on 
the setup (serial or parallel embedding of the dots be- 
tween the source and drain leads), either or both mecha- 
nisms may be significant. In magnetically coupled dots, 
embedded between the leads in series, the conductance is 
low for small exchange coupling when the Kondo singlets 
are formed between each dot and adjacent lead. Con- 
ductance is also low for large exchange coupling, when 
a local singlet state forms between the moments on the 
dots. In contrast, the conductance reaches the unitary 
limiting value of 2e^//i in a narrow interval of J, such 
that J ^ Ta^^'^^- The introduction of additional elec- 
tron hopping between dots breaks the quantum critical 
transition, nevertheless, some signatures of the quantum 
phase transition remain detectable^"^«. 

Strong ferromagnetic RKKY interaction between two 
magnetic impurities coupled to two conduction channels 



leads to three different regimes. At temperatures compa- 
rable to RKKY interaction, ferromagnetic locking of im- 
purity spins occurs; this is followed by a two-stage freez- 
ing out of their local moment as they become screened 
by the conduction electrons^®. This scenario was corrob- 
orated by numerical studies of the two-impurity Kondo 
model^^ and the Alexander- Anderson model2&. Antifer- 
romagnetic and ferromagnetic RKKY interactions lead 
to different transport properties of DQD systems22i^. 
Due to recent advances in nanotechnology, the effects of 
RKKY interaction on transport properties became di- 
rectly observablei^. Conductance through Aharonov- 
Bohm (AB) interferometers with embedded quantum 
dots also depends on the RKKY interactions, which in 
turn depends on the magnetic flus2ii22i2i. A similar 
system of two quantum dots, side-coupled to a single- 
mode channel, allows to study the crossover between fully 
screened and underscreened Kondo impurityM. 



The physics of RKKY interactions is also related to 
the studies of the Kondo effect in integer-spin quantum 
dots^. By tuning the magnetic field, the energy dif- 
ference between singlet and triplet spin states can be 
tuned to zero. At the degeneracy point, a large zero- 
bias resonance with an increased Kondo temperature is 
observed'^^, which can be understood in the framework 
of a two-orbital Anderson model^. 



The interplay of the Kondo effect and the inter- 
impurity exchange interaction leads to a number of inter- 
esting phenomena observed in different realizations of the 
double quantum dot systems. For this reason, we present 
in this work a study of a more general N quantum dot 
systems. Using numerical renormalization group (NRG) 
technique as our primary tool and various analytical ap- 
proaches we investigate the effects of the RKKY inter- 
action in a multi-impurity Anderson model. We present 
results of thermodynamic properties, in particular the 
impurity contribution to the magnetic susceptibility and 
the entropy, as well as various correlation functions. This 
work also provides a setting for further studies of trans- 
port properties of this class of systems. 



The paper is organized as follows. In section |n] we de- 
scribe the class of models under study as well as model 
parameters and approximations used in this work. In 
section IIIII we describe the existence of a hierarchy of 
separated time (and energy) scales and we introduce ef- 
fective models valid at different temperatures. In section 
IIVI we describe the numerical methods that are used in 
section to study the multi- impurity Anderson models. 
Finally, in section IVII we test the stability of the S = 1 
state in the two impurity model with respect to various 
perturbation. Tedious derivations of scaling equations 
and perturbation theory approaches are given in the ap- 
pendices. 



II. THE MODEL 

We study models of N impurities coupled to one single- 
mode conduction channel. The motivation for such mod- 
els comes primarily from experiments performed on sys- 
tems of several quantum dots connected in parallel be- 
tween source and drain electron reservoirs. Since quan- 
tum dots can be made to behave as single magnetic im- 
purities, such systems can be modelled in the first ap- 
proximation as several Anderson impurities embedded 
between two tight-binding lattices as shown schemati- 
cally in Fig. n If the coupling to the left and right 
electrode of each quantum dot is symmetric, it can be 
shown that each dot couples only to the symmetric com- 
bination of conduction electron wave-functions from left 
and right lead, while the antisymmetric combinations 
of wave-functions are totally decoupled and are irrele- 
vant for our purpose^^. We can thus model the parallel 
quantum dots using the following simplified Hamiltonian, 
which we name the "iV-impurity Anderson model" : 

H = Hhand + -ffdots + -ff c • (1) 

X^fco- ^kCka'^ka- is the conduction band 



and 



Here H\^ 

Hamiltonian. Hdots = J2i=i HdotA with 

HdotA 

is the quantum dot Hamiltonian. Finally, 

^c = 4fE(^^4.Cfc. + H.c.) (3) 



SirM~l) + ^{n,^lf ^2) 



Vl 



kai 



is the coupling Hamiltonian, where L is a normaliza- 
tion constant. The number operator rii is defined as 
n.i — ^^di^dicr. Parameter S is related to the more 
conventional on-site energy e^ by (5 = e^ + C//2, where U 
is the on-site Coulomb electron-electron (e-e) repulsion. 
For S — the model is particle-hole symmetric under the 
transformation cj,^ -^ Ck,~(j, dj^ -^ ^di^-a- Parameter & 
thus represents the measure for the departure from the 
particle- hole symmetric point. 

To cast the model into a form that is more convenient 
for a numerical renormalization group study, we make 
two more approximations. We first linearize the dis- 
persion relation ej, of the conduction band, which gives 
Ck — Dk. The wave- number k runs from —1 to 1, there- 
fore 2D is the width of the conduction band. This as- 
sumption is equivalent to adopting a constant density of 
states, po = 1/(213). Second, we approximate the dot- 
band coupling with a constant hybridization strength, 
r = 7rpo|V/cpp. Neither of these approximations affects 
the results in a significant way. In the rest of the pa- 
per, we will present results in terms of the parameters D 
and r, instead of the parameters t and t' of the original 
tight-binding models depicted in Fig. ^ Our notation 
follows that of Refs. I38ll39l for easier comparison of the 
A^-impurity results with the single-impurity case. 
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Figure 1: (Color online) Systems of parallel quantum dots. 
The tight-binding hopping parameter t determines the half- 
width of the conduction band, D = 2i, while parameter t' is 
related to the hybridization F by T/D = {t' /t) . 



where Jk is the effective anti-ferromagnetic Kondo ex- 
change interaction and PqJk = 8T/ttU. This expression 
is vaUd for U <^ D and 5 = 0. 

As we will show later, there is an additional time 
scale Tj ^ ft/JRKKY, originating from the ferromagnetic 
RKKY dot-dot interactions: 



64 r^ 

>/rkky ^ U{pqJk) ~ ^TT7- 



(6) 



From the condition for a well developed Kondo effect, 
?7/(r7r) :^ 1, we obtain Jrkky ^ U . We thus establish 
a hierarchy of time scales tu -^tj <^ tk- 

Based on the three different time-scales, we predict the 
existence of three distinct regimes close to the particle- 
hole symmetric point. The local moment regime is estab- 
lished at T ~ r^*, where T* = U /a and a is a constant of 
the order oneS. In this regime the system behaves as N 
independent spin S — 1/2 impurities. At T ^ Tp, where 
Tp — Jrkky//? and /? is a constant of the order one, 
spins bind into a high-spin S — N/2 state. With further 
lowering of the temperature, at T ~ Tr- the S = N/2 
object experiences the Kondo effect which screens half a 
unit of spin (since there is a single conduction channel) 
to give a ground-state spin of S" — 1/2 = {N — l)/2. 



III. LOW-TEMPERATURE EFFECTIVE 
MODELS 



Schrieffer- Wolff transformation for multiple 
impurities 



Our primary goal is to demonstrate that the low- 
temperature effective model for the multiple impurity 
system is the S = N/2 SU{2) Kondo model: 



H 



/ ^ Jk'k 



Skk' ■ S, 



(4) 



k'k 



where Skk' = hJ2aa'^ka'^aa'Ck'a' is the local-spin den- 
sity in the Wannier orbital in the conduction band that 
couples to all N impurities. S is the collective impu- 
rity S — N/2 spin operator and Jk/k is the momentum- 
dependent anti-ferromagnetic spin-exchange interaction 
that can be derived using the Schrieffer- Wolff transfor- 
mation. Results for J^'k are independent of N. 

We first argue in favor of the validity of the effective 
Hamiltonian, proposed in Eq. Q, by considering the dif- 
ferent time scales of the original A'^-impurity Anderson 
problem. To simplify the argument we further focus on 
the (nearly) symmetric case 5 <^ U within the Kondo 
regime, [//(Ftt) > 1. 

The shortest time scale, tjj ~ h/U, represents charge 
excitations. The longest time scale is associated with 
the Kondo effect (magnetic excitations) and it is given 
by tk ~ h/Tx where Tk is the Kondo temperature of 
the single impurity Anderson models given by Haldane's 
expression 



Tk = 0.182C/v/^;;Jkcxp - 



1 



PqJk 



(5) 



For T < T*, the single impurity Anderson model can 
be mapped using the Schrieffer- Wolff transformation^"^ to 
an s — d exchange model (the Kondo model) with an en- 
ergy dependent anti-ferromagnetic exchange interaction 
Jk'k- In this subsection we show that for multiple impuri- 
ties a generalized Schrieffer- Wolff transformation can be 
performed and that below T* , the A^-impurity Anderson 
model maps to the A^-impurity S* = 1/2 Kondo model. 
Furthermore, the exchange constant is shown to be the 
same as in the single impurity case. 

Due to the hybridization term Vfe, the electrons are 
hopping on and off the impurities. Since all impuri- 
ties are coupled to the same Wannier orbital, it could 
be expected that these hopping transitions would some- 
how "interfere" . It should be recalled, however, that the 
dwelling time tu is much shorter than the magnetic time 
scales Tj and tk- In other words, spin- flips are realized 
on a much shorter time-scale compared to the mean-time 
between successive spin-flips; for this reason, each local 
moment may be considered as independent. Note that 
the impurities do in fact "interfere" : there are 0{V^) pro- 
cesses which lead to an effective ferromagnetic RKKY ex- 
change interaction between pairs of spins and ultimately 
to the ferromagnetic ordering of spins at temperatures 
below '^ Jrkky- This will be discussed in the following 
subsection. 

The Schrieffer- Wolff transformation is a canonical 
transformation that eliminates hybridization terms Vk to 
first order from the Hamiltonian H, i.e. it requires that^ 



H = e^He- 



(7) 



have no terms which are first order in Vk- We expand H 
in terms of nested commutators: 



H = H+[S,H] + ^[S,[S,H]] 



(8) 



and write H = Ho + He, where Hq = i^band + ^dots- We 
then choose S to be first order in Vk so that 



[S, Ho] +H, = 0. 



(9) 



As previously discussed, each impurity can be considered 
independent due to the separation of time scales. There- 
fore, we choose the generator S to be the sum S = J2i ^i 
of generators Si , where the generator St for each impurity 
has the same form as in the single-impurity case: 






kaa 



Efe — Eq 



H,-a'^kcr"'^cr 



H.c. 



(10) 



with e± = 6zLU/2 and the projection operators nf_^ are 
defined by 



n7 „ = 1 - rii _„ 



(11) 



The resulting effective Hamiltonian is then given by 

H,s^Ho + ^[S,H,], (12) 

which features 0{VJ^) effective interactions with the lead- 
ing terms that can be cast in the form of the Kondo an- 
tiferromagnetic exchange interaction 



i \ kk' 



■.^kk' ' ^i \ 1 



(13) 



where Si is the S — 1/2 spin operator on impurity i 
defined by 5'^ = ^ X]qq' di^crdia' and the exchange con- 
stant Jk'k is given by 



Jk'k=Vk,Vk{ + 



1 



ek-i6 + U/2) ek' -{5 + U/2) 



ek-{5- U/2) ek' -{5- C//2) 



(14) 



If we limit the wave- vectors to the Fermi surface, i.e. for 
k ~ k' ^ kp, 'WG obtain 



Jk = 2\Vk, 



1 



1 



\S-U/2\ \S + U/2\ 



(15) 



This result is identical to Jk'k obtained for a single 
impurity^. 

As it turns out, the Schrieffer- Wolff transformation, 
Eqs. lf7|l- H12|) . produces inter-impurity interaction terms 



in addition to the expected impurity-band interaction 
terms. In the particle- hole symmetric case (i5 = 0), these 
additional terms can be written as 



AH, 




where 



cff 



i<j,cr 



N h 



hop J 



d]„d.ic 



(16) 



(17) 



Since the on-site charge repulsion favors states with single 
occupancy of each impurity, the term in the parenthesis 
in Eq. (|16|) is on the average equal to zero. Furthermore, 
if each site is singly occupied, possessing small fluctua- 



tions of the charge {nf 



0, hopping between the 



sites is suppressed and the term /ihop represents another 
small factor. The Hamiltonian AH is thus clearly not 
relevant: impurities are indeed independent. 

On departure from the particle-hole symmetric point 
((5^0), Ai/off generalizes to 



Ai/eff = 2 



U\Vk? 

C/2 _ 4^2 



Y^n,^N]-2N-]h 



hop- 



(18) 
For moderately large S/U this Hamiltonian term still rep- 
resents only a small correction to Eq. IjlSI) . However, 
for strong departure from the particle-hole (p-h) sym- 
metric point, close to the valence-fluctuation regime (i.e. 
5 -^ U/2), the A_ffcff becomes comparable in magnitude 
to Hex and generates hopping of electrons between the 
impurities. 

The above discussion leads us to the conclusion that 
just below T* the effective Hamiltonian close to the p-h 
symmetric point is 



H 



cff 



— Hhund + 2_^ 2^ Jk'kSkk' ' Si 



(19) 



k'k 



If the dots are described by unequal Hamiltonians H,iot,i 
or have unequal hybridizations V/:, then the mapping of 
the multi-impurity Anderson model to a multi-impurity 
Kondo model still holds, however with different effective 
exchange constants J^,^.. 



B. RKKY interaction and ferromagnetic spin 
ordering 

We now show that the effective RKKY exchange in- 
teraction between the spins in the effective A^-impurity 
Kondo model, Eq. H19|l . is ferromagnetic and also respon- 
sible for locking of spins in a state of high total spin for 
temperatures below T < Jrkky- 

The ferromagnetic character of the RKKY interaction 
is expected, as shown by the following qualitative argu- 
ment. We factor out the spin operators in the effective 



Hamiltonian Eq. (fT^ : 



^cff = -H^band + I 2_^ Jk'k^kk' 1 ' 2^ Si 



(20) 



,k'k 



Spins Si are aligned in the ground state since such orien- 
tation minimizes the energy of the system. This follows 
from considering a spin chain with N sites in a "static 
magnetic field" X^fe'fc Jk'k^kk' ■ The assumption of a static 
magnetic field is valid due to the separation of relevant 
time scales, tk ^ tj. States with S < N/2 are clearly 
excited states with one or several "misaligned" spins. 

Since the inter-dot spin-spin coupling is a special case 
of the Ruderman-Kittel-Kasuya-Yosida (RKKY) interac- 
tion in bulk systems^^ , a characteristic functional depen- 
dence given by 



-'^RKKY ocU{poJKf 



64 r^ 

^77 



U 



(21) 



is expected. The factor U in front of {pqJkY plays 
the role of a high-energy cut-off, much like the 0.196f/ 
effective-bandwidth factor in the expression for Tk^ 
Eq. 0. 

Using the Rayleigh-Schrodinger perturbation theory 
we calculated the singlet and triplet ground state ener- 
gies Es and Et to the fourth order in Vk for the two- 
impurity case (see Appendix 0. We define the RKKY 
exchange parameter by Jrkky — Es~Et] positive value 
of Jrkky corresponds to ferromagnetic RKKY interac- 
tion. For U/d < 0.1, the prefactor of (poJk)^ in the 
expression (I21II is indeed found to be linear in U. To- 
gether with the prefactor the perturbation theory leads 
to 



Jrkky = 0.62U{poJk)^ for U/D « 1, 



(22) 



which, as we will show later, fits very well our numerical 
results. The RKKY interaction becomes fully established 
for temperatures below Tj which is roughly one or two 
orders of magnitude smaller than T* (Tj is defined in 
Appendix 0. Since the RKKY interactions in the first 
approximation do not depend much on the number of im- 
purities, for A^ > 2 the exchange interaction between each 
pair of impurities has the same strength as in the two im- 
purity case. Therefore, for temperatures just below Tj, 
the effective Hamiltonian for the A^-impurity Anderson 
model becomes 



In conclusion, for T <Tp the states with total spin less 
than A^/2 can be neglected, and the system behaves as 
if it consisted of a single spin S of magnitude S — N/2. 
The effective Hamiltonian at very low temperatures is 
therefore the S = A/2 SU{2) Kondo model 



H, 



cff 



E^- 



kSk,k' 



(24) 



k.k' 



where S = V {J2i Si)V and V is the projection operator 
on the subspace with total spin S — N/2. Other multi- 
plets are irrelevant at temperatures below Tp. We point 
out that the Kondo temperature for this model is given 
by the formula for the single impurity Anderson model, 
Eq. 0, irrespective of the number of dots A, since the 
ferromagnetic interaction only leads to moment order- 
ing, while the exchange interaction of the collective spin 
is still given by the same Jk'k- 

It should be mentioned that if the exchange constants 
•^k'k ^'^^ different impurities are different, there will be 
some mixing between the spin multiplets. The simple 
description of impurities as a collective S = N/2 spin still 
holds even for relatively large differences, but in general 
the virtual excitations to other spin multiplets must be 
taken into account. This is studied in detail for the case 
of two dots in Section FV I PI 



IV. THE METHODS 

A. Numerical renormalization group 

The method of choice to study the low-temperature 
properties of quantum impurity models is the Wilson's 
numerical renormalization group (NRG)22i22i^. The 
NRG technique consists of logarithmic discretization of 
the conduction band described by iJband, mapping onto 
a one-dimensional chain with exponentially decreasing 
hopping constants, and iterative diagonalization of the 
resulting Hamiltonian. Since all A^ impurities couple to 
the band in the same manner, they all couple to the same, 
zero-th site of the chain Hamiltonian^^: 



2.^2-^ ^" /n,cr/n+l,o- + /n+l,o-/n 



n=0 o- 



(25) 



ff, 



cff 



■^d+ [Y. Jk'kSkkA E •^-•^RKKY E ^--^y + iJdot« + E (^) (/ol'^- + Aahc 



(23) 
When the temperature drops below a certain tempera- 
ture Tp, the spins align and form a ferromagnetically- 
frozen state of maximum spin S — N/2. The transition 
temperature Tp is generally of the same order as Jrkky, 
i.e. Tp = Jrkky//?, where /3 is an A^-dependent constant 
of the order one. This relation holds if Tp <C Tj, other- 
wise Tp needs to be determined using a self-consistency 
equation ljA5|l . as discussed in Appendix \K\ 



Here /^^ are the chain creation operators and ^„ are 
constants of order 1. In addition to the conventional 
Wilson's discretization scheme^, we also used Campo 
and Oliveira's new discretization approach using an over- 
complete basis of states^ with A = 4, which improved 
convergence to the continuum limit. We made use of the 
"2;-trick" with typically 6 equally spaced values of the 
parameter d^. 



1. Symmetries 

The Hamiltonian Q has the foUowing symmetries: 
a) C/(l)gauge symmetry due to global phase (gauge) in- 
variance. The corresponding conserved quantity is the 
total charge (defined with respect to half-filling case): 
Q = ^j (rii — 1), where the sum runs over all the impu- 
rity as well as the lead sites; b) SU(2)spin spin symmetry 
with generators S = J2i^Y.aa' ^la^^aa'aia', where cr 
are the Pauli matrices. Since operators Q, S^ and 5*2 
commute with H, the invariant subspaces can be classi- 
fied according to quantum numbers Q, S and Sz- Com- 
putation of matrix elements can be further simplified us- 
ing the Wigner-Eckart theorem^. 

In the particle-hole symmetric point, i.e. S = 0, Hamil- 
tonian has an additional SU{2)iso isospin symmetryii. 
We define isospin operators on impurity site i using 



J2 ^'L' 



(26) 



J 



where the Nambu spinor rjl on the impurity orbitals is 
defined by 



In all our NRG calculations we took into account the 
conservation of the charge and the rotational invariance 
in the spin space, i.e. the f7(l)gaugc (^ SU{2)spin symme- 
try which holds for all perturbed models considered, or 
the 5C/(2)spin(8iS'[/(2)isospin symmetry where applicable. 
The number of states that we kept in each stage of the 
NRG iteration depended on the number of the dots N, 
since the degeneracy increases exponentially with N: ap- 
proximately as 4^ at the high-temperature free orbital 
regime and as 2^ in the local-moment regime. In the 
most demanding A^ = 4 calculation we kept up to 12000 
states at each iteration (which corresponds to > 32000 
states taking into account the spin multiplicity of states), 
which gave fully converged results for the magnetic sus- 
ceptibility. 

For large scale NRG calculations it is worth taking 
into account that the calculation of eigenvalues scales 
as 0{n^) and the calculation of eigenvectors as 0{n^), 
where n is the dimension of the matrix being diagonal- 
ized. Since eigenvectors of the states that are truncated 
are not required to recalculate various matrices prior to 
performing a new iteration, considerable amount of time 
can be saved by not calculating them at all. 



Vi 



= ^^A 



(27) 



We also define 1+ = r + iiv and /- = (/+)t. We 
then have, for example, I[ — [rii — l)/2 = Qi/2 and 
I^ = d\^d\^. The isospin symmetry is thus related to the 
electron pairing. In terms of the isospin operators the 
impurity Hamiltonian can be written as 



H, 



dot.1 



2«f + AU{I, 



z\2 



25If 



\ui,.f. 



(28) 



where we took into account that for spin- 1/2 operators 
(Pauh matrices) (/f)^ = 1/3(7^)2. 

On the Wilson chain the isospin is defined similarly 
but with a sign alternation in the definition of the Nambu 
spinors ^„: 



a 



■In,] 
(-l)"/n,i 



(29) 



The total isospin operator is obtained through a sum 
of li for all orbitals of the problem (impurities and con- 
duction band). For (5 = 0, both I and I^ commute with 
H and / and I^ are additional good quantum numbers. 
Note that /^ = Q/2, therefore C/(l)gaugc is in fact a sub- 
group of SU {2.)iso- Due to isotropy in isospin space, the 
Iz dependence can again be taken into account using the 
Wignert-Eckart theorem. 

Spin and isospin operators commute, [Si,Ij\ = 
for all i,j. Therefore, for (5 = the problem has a 
SU(2)s-pin ® SU{2)[so symmetry which, when explicitly 
taken into account, leads to a further significant reduc- 
tion of the numerical task. 



2. Calculated quantities 

We have computed the following thermodynamic quan- 
tities 

• the temperature-dependent impurity contribution 
to the magnetic susceptibihty XimpiT) 



Ximp(r) 



knT 



{{S'z) - {S'zk 



(30) 



where the subscript refers to the situation when 
no impurities are present (i.e. H is simply the band 
Hamiltonian -ffband), 9 is the electronic g factor, 
/is the Bohr magneton and ks the Boltzmann's 
constant. It should be noted that the combination 
^Ximp/(5M-B)^ can be considered as an effective mo- 
ment of the impurities, ^eff- 

• the temperature-dependent impurity contribution 
to the entropy S-mipiT) 



'-'impi-' J 



(E-F) (E-F), 



T 



T 



(31) 



where E = (H) = Ti {Her"'^^^^^) and 
F = -fcBTlnTr(e"^/''B^). From the quantity 
Simp/kB we can (deduce the effective degrees-of- 
freedom v of the impurity as 5'i,np/fcB ~ Inz^. 

• thermodynamic expectation values of various op- 
erators such as the on-site occupancy (rij), local 
charge-fiuctuations (((5n)^) = (nf) — (rti)^, local- 
spin (S'j) and spin-spin correlations {Si ■ Sj). 



In the following we drop the sufRx imp in Ximp, but 
one should keep in mind that impurity contribution to 
the quantity is always implied. We also set A;b = 1- 



B. Bethe Ansatz 

The single-channel SU{2) Kondo model can be ex- 
actly solved for an arbitrary spin of the impurity using 
the Bethe Ansatz (BA) methodMi^^i^. This technique 
gives exact results for thermodynamic quantities, such 
as magnetic susceptibility, entropy and heat capacity. It 
is, however, incapable of providing spectral and trans- 
port properties. For the purpose of comparing results 
of the single-channel SU{2) Kondo model with NRG re- 
sults of the N-impurity Anderson model, we have nu- 
merically solved the system of coupled integral equations 
using a discretization scheme as described, for example, 
inRef.Ei 



C. Scaling analysis 

Certain aspects of the Kondo physics can be correctly 
captured using the perturbative renormalization group 
approach based on the "poor-man's scaling" technique 
due to Anderson^^. A brief account of this method is 
given in Appendix IbI 



V. NUMERICAL RESULTS 

We choose the parameters U and F so that the rele- 
vant energy scales are well separated which enables clear 
identification of various regimes and facilitates analytical 
predictions. 

In Fig. [5] we show temperature dependence of magnetic 
susceptibility and entropy for N = 1,2,3 and 4 systems. 
As the temperature is reduced, the system goes through 
the following regimes: 

1. At high temperatures, T > Tj*, the impurities are 
independent and they are in the free orbital regime 
(FO) (states |0), | t), | i) and |2) on each impurity 
are equiprobable) . Each dot then contributes 1/8 
to HeS = Tx/igfJ-s)^ for a total of HeS = N/8. The 
entropy approaches S'imp = iV In 4 since all possible 
states are equally probable^. 

2. For Tp < T < T* each dot is in the local-moment 
regime (LM) (states | t) and | |) are equiprobable, 
while the states |0) and |2) are suppressed). Each 
dot then contributes 1/4 to fieS for a total of iV/4. 



The entropy decreases to Sh 



TV In 2. 



3. For Tk < T < T^ and N > 1 the dots lock into 
a high spin state S = N/2 due to ferromagnetic 
RKKY coupling between local moments formed 
on the impurities. This is the ferromagnetically 



frozen regime (FF)2& with /XeS = S{S + l)/3 = 
N/2{N/2 + l)/3. The entropy decreases further to 



Oil 



ln(2S'+l) =ln(iV+l). 



4. Finally, for T < Tk, the total spin is screened from 
S = N/2 to S ^ S -1/2 = {N - l)/2 as we en- 
ter the partially-quenched, Kondo screened strong- 
coupling (SC) A^-impurity regime with jj^eff = S{S-\- 
l)/3 = (TV - l)/2[(iV - l)/2 + l]/3. The remaining 
S — 1/2 spin is a complicated object: a S = N/2 
multiplet combination of the impurity spins anti- 
ferromagnetically coupled by a spin-1/2 cloud of 
the lead^^. In this regime, the entropy reaches its 



minimum value of S'imp — ln(2iS' + 1) = In TV. 
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Figure 2: (Color online) a) Temperature-dependent suscepti- 
bility and b) entropy of the A*'-dot systems calculated using 
the NRG. The symbols in the susceptibility plots were cal- 
culated using the thermodynamic Bethe Ansatz approach for 
the corresponding 5" = N/2 SU(2) Kondo models (• S = 1/2, 
■ 5 = 1, ♦ 5 = 3/2, A S = 2). 

In Fig. 121 atop the NRG results we additionally plot 
the results for the magnetic susceptibility of the S — N/2 
SU{2) Kondo model obtained using an exact thermody- 
namic Bethe-Ansatz method. For T < Tp nearly per- 
fect agreement between the A'^-impurity Anderson model 
and the corresponding S = N/2 SU{2) Kondo model are 
found over many orders of magnitude. ^^ This agreement 



N_ 

1 

2 

3 

4 



Kondo temperature Tk/D LM-FO temperature T'^/D ^he transition from to 1/4 is realized^ at T 
77^ ca,n extra,ct Tf, iisinff the (somewhat a,rbitra,rv 



1.20 X 10" 
1.23 X 10" 
1.29 X 10" 
1.32 X 10" 



T^. We 
can extract Tp using the (somewhat arbitrary) condition 



1.87 X lO"'' 
2.11 X lO"'^ 
2.32 X 10"^ 



Table I: Kondo temperatures for different numbers of quan- 
tum dots A'^ corresponding to plots in Fig. |5| 



is used to extract the Kondo temperature of the multiple- 
impurity Anderson model. The fitting is performed nu- 
merically by the method of least-squares; in this manner 
very high accuracy of the extracted Kondo temperature 
can be achieved. The results in TableQlpoint out the im- 
portant result of this work that the Kondo temperature 
is nearly independent of iV, as predicted in Section IlIIBI 
In this sense, the locking of spins into a high-spin state 
does not, by itself, weaken the Kondo effecl!i2i2£; how- 
ever, it does modify the temperature-dependence of the 
thermodynamic and transport properties^^*^. 

It is instructive to follow transitions from high- 
temperature FO regime to LM and FF regime through a 
plot combining the temperature dependence of the mag- 
netic susceptibility and of other thermodynamic quanti- 
ties, as presented in Fig.|21 Charge fluctuations {{Sn)^) 
show a sudden drop at T ^ Tj* representing the FO - LM 
transition. In contrast, the magnitude of the total spin 
S increases in steps: 5* = 1/2, (-s/T — l)/2 and 1. Val- 
ues of S in these plateaus are the characteristic values of 
doubly occupied double-quantum dot system in the FO, 
LM and FF regime, respectively. 
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Figure 3: (Color online) Temperature-dependence of suscep- 
tibility, charge fluctuations {{5n)^), total spin S and the spin- 
spin correlations {Si ■ S2) of the 2-dot system. 

The LM-FF transition temperature Tp can be deduced 
from the temperature dependence of the spin-spin corre- 
lation function. In the FF regime the spins tend to align, 
which leads to (Si • S2) ^- 1/4 as T ^ 0, see Fig. El 



(Si-S2)(r*) = i/2(Si.S2)(r^o). 



(32) 



In section IVIB II we show that this condition is in very 
good agreement with Tp — Jbxky I P obtained by deter- 
mining the explicit inter-impurity antiferromagnetic cou- 
pling constant J12, defined by the relation Jrkky + >/i2 = 
that destabilizes the high-spin S = N/2 state. The 
extracted Tp transition temperatures that correspond 
to plots in Fig. |21 are given in Table |l| We find that 
they weakly depend on the number of impurities, more 
so than the Kondo temperature. The increase of Tp 

for 
for N 



with N can be partially explained by calculating Tp 



a spin Hamiltonian H = — Jrkky X]i<i "^j ' ^j 

spins decoupled from leads. Using Eq. H32() we obtain 

T* 



1.18 Jrkky for iV = 2, T* 



1.36 Jrkky for 



TV = 3 and r|^ w 1.55 Jrkky for TV = 4. 

By performing NRG calculations of Tp for other pa- 
rameters U and r and comparing them to the predic- 
tion of the perturbation theory, we found that the simple 
formula (|22|) for Jrkky agrees very well with numerical 
results. 

The effect on thermodynamic properties of varying U 
while keeping T/U (i.e. pqJk) fixed is illustrated in Fig.^ 
for 2- and 3-dot systems. Parameters F and U enter ex- 
pressions for Tp = Jrkky//? and Tk only through the ra- 
tio T/U, apart from the change of the effective bandwidth 
proportional to U, see Eq. and H22I) . This explains 
the horizontal shift towards higher temperatures of sus- 
ceptibility curves with increasing U, as seen in Fig. 2Jl. 
The NRG results and the Bethe-Ansatz for the Kondo 
models with S = 1 and 5' = 3/2 show excellent agree- 
ment for T < Tp. In Figs.^JD and^l: we demonstrate the 
nearly linear [/-dependence of Tp and T^, respectively. 

In Fig. Owe show the effect of varying T/U while keep- 
ing U fixed. In this case, T* stays the same, Tp is shifted 
quadratically and Tk exponentially with increasing T/U . 
Fig. |S1d shows the agreement of Tp with expression H22|l , 
while Fig. \^ shows the agreement of the extracted values 
of Tk with formula (O . 

We note that for A^ > 2, eventual coupling to an ad- 
ditional conduction channel (for example, due to a small 
asymmetry in the coupling to the source and drain elec- 
trodes) would lead to screening by additional half a unit 
of spin^^-^^ and the residual ground state spin would be 
S -1 = N/2 - 1. For TV > 3 and three channels (due to 
weak coupling to some third electrode), three half-units 
of spin would be screened, and so forth. These addi- 
tional stages of Kondo screening would, however, occur 
at much lower temperatures; all our findings still apply 
at temperatures above subsequent Kondo cross-overs. 

In systems of multiple quantum dots, an additional 
screening mechanism is possible when after the first 
Kondo cross-over, the residual interaction between the 
remaining spin and the Fermi liquid quasi-particles is 
antiferromagnetic^^. This leads to an additional Kondo 
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Figure 4: (Color online) a) Temperature-dependent suscepti- 
bility of the 2 and 3-dot systems with the same T /U ratio. 
Open (filled) symbols are Bethe-Ansatz results for the S = 1 
(S = 3/2) Kondo model, b) Comparison of LM-FF transition 
temperature Tp with predictions of the perturbation theory, 
c) Comparison of calculated Tk with the Haldane's formula. 
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Figure 5: (Color online) a) Temperature-dependent suscepti- 
bility of the 2-dot system for equal e-e repulsion U jD = 0.01 
and for different hybridization strengths P. Symbols repre- 
sent the Bethe-Ansatz susceptibility for the S = \ Kondo 
model with corresponding Tk- b) Comparison of calculated 
and predicted Tp. c) Fit of Tk to the Haldane's formula, 
Eq. 0. 



cross-over at temperatures that are exponentially smaller 
than the first Kondo temperature. Such two-stage Kondo 
effect occurs, for example, in side-coupled double quan- 
tum dots^fliSiiSiS and triple quantum dot coupled in 
series^. In parallelly coupled systems, the residual inter- 
action between the remaining spin and the Fermi liquid 
quasi-particles is, however, ferromagnetic as can be de- 
duced from the splitting of the NRG energy levels in the 
strong-coupling fixed point "^^i the strong-coupling fixed 
point is stable. 

We have thus demonstrated that with decreasing tem- 
perature the symmetric (b = 0) multi-impurity Ander- 
son model flows from the FO regime, through LM and 
FF regimes, to a stable underscreened S = N jl Kondo 
model strong-coupling fixed point. The summary of dif- 
ferent regimes is given in Table ITU 



VI. STABILITY OF TV = 2 SYSTEMS WITH 
RESPECT TO VARIOUS PERTURBATIONS 

We next explore the effect of various physically rele- 
vant perturbations with a special emphasis on the ro- 
bustness of the ferromagnetically frozen state and the 
ensuing S — N/2 Kondo effect against perturbation of 
increasing strength. We show that the system of multi- 
ple quantum dots remains in a 5* = N /2 state even for 
relatively large perturbations. We also study the quan- 
tum phase transitions from the S = N/2 state driven 
by strong perturbations. In this sections we limit our 
calculations to the N = 2 system. 



A. Variation of the on-site energy levels 

1. Deviation from the particle-hole symmetric point 

A small departure from the particle-hole symmetric 
point ((5 7^ 0) does not destabilize the S = N/2 Kondo 
behavior: the magnetic susceptibility curves still follow 
the Bethe-Ansatz results even for S/U as large as 0.4, 
see Fig. |S^. For 5 > Sc, where Sc/D ^ 0.45 is the 
critical value of parameter (5, the triplet state is desta- 
bilized. Consequently, there is no Kondo efect. This 
is a particular case of the singlet-triplet transition that 
is a subject of intense studies in recent years, both 
experimentally^^i^^i^''" and theoreticallji^SiSiSiSiffi. 

In the asymmetric single impurity model, the valence- 
fiuctuation (VF) regime is characterized by HeS = 
Tx{T)/igf^B)^ ~ 1/622.. The VF regimes occurs at T* 
and the transition from VF to LM regime occurs at 
T2* ^ \E^\/a, where E^ is the renormalized on-site en- 
ergy of the impurity: E^ ^ ej,- ^\n {-U/E^). For two 
uncorrelated dots in the VF regime, we expect /icff ^1/3. 
In Fig. |HJi we plotted a number of susceptibility curves 
for parameters S in the proximity of the singlet-triplet 
transition. While there is no clearly-observable valence- 
fiuctuation plateau, the value of ficS is indeed near 1/3 



Table II: Regimes of the symmetric (5 = 0) A''-impurity Anderson model 
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Regime 


Relevant states 


Magnetic susceptibility 


Spin correlations 


Charge fluctuations 


Entropy 






Moff = TXnr.p{T)/{gilBf 


{Si -Si) 


in') - {nf 


•^imp 


FO 


iVx(|0),|T),U>,|2)) 


N/8 





0(1) 


NlnA 


LM 


^x(|T),U» 


N/4: 





small 


Nln2 


FF 


\S = N/2,S,) 


N/2{N/2 + l)/3 


-1/4 


small 


ln{N + 1) 


SC 


\S = N/2-l/2,S,) 


(iV - l)/2(iV/2 + l/2)/3 


-1/4 


small 


\nN 




0.182|£^^|. The Kondo temperature is now given by 



10 0.1 0.2 0.3 0.4 0.5 
5/U 



% 0.1 0.2 0.3 0.4 0.5 



5/U 



Figure 6: (Color online) a) Temperature-dependent suscep- 
tibility of the 2-dot systems on departure (5 > 0) from the 
particle-hole symmetric point, 5 = 0. Symbols are fits to 
the universal susceptibility obtained using the Bethe-Ansatz 
method for the S = 1 Kondo model, b) Calculated and pre- 
dicted Kondo temperature, Eq. 1351 . For comparison we also 
plot Tk given by Eq. 0, which shows expected discrepancy 
for large S/U. c) Calculated Tp and the fit to an exponential 
function. 



between T^ and Ti{Sc). 

In Fig. 03 we compare calculated Kondo temperatures 
with analytical predictions based on the results for the 
single impurity model'^^. For (5^0, Jk generalizes ac- 
cording to the SchriefFer- Wolff transformation, Eq. (|15|l . 
Departure from the p-h symmetric point also induces po- 
tential scattering 



PqK 



1 



T 



2tt \\5-U/2\ \5 + U/2\ 



(33) 



The effective Jk that enters the expression for the Kondo 
temperature i^ 



Jk — Jk 



1 + [ttp^KY 



(34) 



Tk = 0.182|j;^|VpoJKexp -1/(poJk) 



(35) 



This analytical estimate agrees perfectly with the NRG 
results: for moderate 5/U , the results obtained for asym- 
metric single impurity model also apply to the multi- 
impurity Anderson model. 

In Fig. IHt we show the (5-dependence of the LM-FF 
transition temperature Tp. Its value remains nearly in- 
dependent of 5 in the interval 6 < 0.4C/ and then it sud- 
denly drops. More quantitatively, the dependence on 5 
can be adequately described using an exponential func- 
tion 



T*p{5) = t;(o) 



1 — exp 



5~5c 
A 



(36) 



where Tp{0)/D = 1.8 x 10~^ is the transition temper- 
ature in the symmetric case, Sc/D = 0.45 is the criti- 
cal S and X/D = 2.1 x 10^^ is the width of the transi- 
tion region. Exchange interaction Jrkky does not de- 
pend on S for U/D = 0.01 ^ 1, which explains con- 
stant value of Tp{S) for S < QAU. At a critical value 
Sc, Tp goes to zero and for still higher S the spin- 
spin correlation becomes antiferromagnetic. Since the 
ground-state spins are different, the triplet and singlet 
regime are separated by a quantum phase transition at 
S = Sc- This transition is induced by charge fluctuations 
which destroy the ferromagnetic order of spins as the sys- 
tem enters the VF regime. The exponential dependence 
arises from the grand-canonical statistical weight factor 
exp[(5(n — 2)/{kBT)], where n is the number of the elec- 
trons confined on the dots. The transition is of the first 
order, since for equal coupling of both impurities to the 
band there is no mixing between the n — 2 triplet states 
and the n = singlet stated. 

For S slightly lower than the critical Sc, the effective 
moment Tx{T) shows a rather unusual temperature de- 
pendence. It first starts decreasing due to charge fiuc- 
tuations, however with further lowering of the temper- 
ature the moment ordering wins over, Tx{T) increases 
and at low-temperatures approaches the value charac- 
teristic for the partially screened S — 1 moment, i.e. 



and the effective bandwidth 0.182C/ is replaced by Tx/igfJ-s) '^ 1/4- 



11 



2. Splitting of the on-site energy levels 

We next consider the 2-dot Hamiltonian with unequal 
on-site energies Sf. 



Hdot,i = 6i{nt - 1) + 77(^1 - 1) 



(37) 



We focus on the case 61 = A and 62 = —A, which 
represents another experimentally relevant perturbation. 
This model is namely particle-hole symmetric for an arbi- 
trary choice of A under a generalized p-h transformation 



Cfc. 



dl 



d-2 



A. 



dl, 



The total oc- 



cupancy of both dots is exactly 2 for any A. We can 
therefore study the effect of the on-site energy splitting 
while maintaining the particle-hole symmetry. Suscepti- 
bility curves are shown in Fig. [7^ for a range of values 
of A. For A up to some critical value Ac ~ 0.47 the 
2-dot Anderson model remains equivalent to the S — 1 
Kondo model for T < Tp. A singlet-triplet transition of 
the Kosterlitz-Thouless typo^°i^^ occurs at A = Ac. 



Due to the particle-hole symmetry no potential scattering 
is generated. The effective Kondo Hamiltonian for small 
A is thus nearly the same as that for small S discussed in 
the previous section. In Fig. [7|d calculated Kondo tem- 
peratures are plotted in comparison with analytical result 
from Eqs. (0) and H38(l . The agreement is excellent. 

The properties of the systems with (5^0 and A 7^ 
become markedly different near respective singlet-triplet 
transition points. For 5^0, the transition is induced 
by charge fluctuations which suppress magnetic order- 
ing and, due to equal coupling of both dots to the band, 
the transition is of first order. For A y^ the transition 
is induced by depopulating dot 2 and populating dot 1 
while the total charge on the dots is maintained, which 
leads to the transition from an inter-impurity triplet to a 
local spin-singlet on the dot 1 . Since there is an asymme- 
try between the dots, the transition is of the Kosterlitz- 
Thouless type^. 

The Kondo temperature of the S — 1/2 Kondo screen- 
ing near the transition on the singlet side, T*, is approx- 
imately given by 
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Figure 7: (Color online) a) Temperature-dependent suscep- 
tibility of the 2-dot system with unequal (detuned) on-site 
energies, Si — A, S2 = —A. Full symbols present Bethe- 
Ansatz results of the equivalent S = 1 Kondo model, while 
empty symbols are BA results of a S = 1/2 Kondo model, b) 
Comparison of calculated and predicted Kondo temperature, 
see Eqs. @ and 13811 . c) The Kondo temperature of the 
S=l/2 Kondo screening on the singlet side of the transition 
and a fit to Eq. 1391 . 

Even though the two dots are now inequivalent, the 
Schrieffer- Wolff transformation yields the same Jk for 
both spin impurities. We obtain 



Jk = 2\Vk, 



1 



1 



U/2\ 



U/2\ 



(38) 



logT7i:> = -a-/3exp 



A-A 



(39) 



We obtain a « 7, /3 « 2.8, A/D « 0.477 and X/D « 1.5 x 
10"'^. This expression is consistent with the cross-over 
scale formula T* oc cxp[— rR-/Ji2] for a system of two 
fictitious spins, one directly coupled to the conduction 
band and the other side-coupled to the first one with 
exchange-interaction J12 that depends exponentially on 
A: Ji2=TK//3exp[(A-A)/A]. 



B. Inter-impurity interaction 

1. Inter-impurity exchange interaction 

In this subsection we show that by introducing an 
explicit exchange interaction J12 between the localized 
spins on the dots, the strength of the RKKY interaction, 
^RKKY, can be directly determined. We thus study the 
two-impurity Anderson model with 



H. 



dots 



-E^<^ 



ot.i 



J12S1 ■ S2, 



where J12 > 0. 

As seen from Fig. |H1 for J12 above a critical value Jc, 
the RKKY interaction is compensated, local moments on 
the dots form the singlet rather than the triplet which in 
turn prevents formation of the S — 1 Kondo effect. The 
phase transition is of the first order^. Using Eq. (|32|) , we 
obtain Tp/D « 1.87 x 10^^ for the non-perturbed prob- 
lem with the same U and F, while Jc/D « 1.68 x 10~^. 
Taking into account the definition Tp = Jrkky / f3, where 
/3 ~ 1 , we conclude that Jrkky agrees well with the crit- 
ical value of Jc, i.e. Jc — Jrkky- The perturbation 
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theory prediction of Jkkky /D — 1.6 x 10 ^ also agrees 
favorably with numerical results. 
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Figure 8: (Color online) Temperature-dependent susceptibil- 
ity of the 2-dot systems for different anti-ferromagnetic inter- 
impurity couplings Ji2. Circles are Bethe-Ansatz results for 
the susceptibility of the S = \ Kondo model with the Kondo 
temperature which is equal for all cases where J12 < Jc. 

As long as J12 < Jc, even for J12 > Tr-, the S = \ 
Kondo effect survives and, moreover, the Kondo temper- 
ature remains unchanged, determined only by the value 
of po^K as in the J\i — case. The only effect of increas- 
ing J12 in the regime where J\i < Jc is the reduction of 
the transition temperature into the triplet state, which 
is now given by Tp ~ Jcs/ P with the effective inter- 
impurity interaction J^q = Jrkky — Ji2- 



2. Hopping between the impurities 

We now study the two-impurity Anderson model with 
additional hopping between the dots: 



H, 



dots 



E 



Hdot,i - tl2 ^ (4o-^2<T + 4o-^l 



(40) 



This model can be viewed also as a single-channel ver- 
sion of the Alexander- Anderson modelSi in the limit of 
zero separation between the impurities. The magnetic- 
susceptibility curves are shown in Fig. ^ 

The hopping leads to hybridization between the atomic 
levels of the dots which in turn results in the formation 
of an even and odd level ( "molecular orbital" ) with ener- 
gies ee,o = Ed i ii2. In the presence of interaction U 
there are two contributions to the energy of the low- 
lying states: "orbital energy" proportional to ii2 and 
"magnetic energy" due to an effective antiferromagnctic 
exchange Jafm = 4t^2/t^j which is second-order in ii2. 
Even though the orbital energy is the larger energy scale, 
the Kondo effect is largely insensitive to the resulting 
level splitting. Instead, the Kondo effect is destroyed 
when Jafm exceeds Jrkky : much like in the case of ex- 
plicit exchange interaction between the dots which was 
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Figure 9: (Color online) Temperature-dependent susceptibil- 
ity of the 2-dot systems with inter-dot tunneling coupling ii2 . 
For ti2/D « 2 X 10"", we have Jafm/D « 1.6 x 10"^ which 
agrees well with the critical value of Jc/D ~ 1.7 x 10~^ found 
in the case of an explicit exchange interaction between the 
dots, see Fig.|H| 



discussed in the previous subsection. We should empha- 
size the similarity between the curves in Figs. |H1 and El 

In the wide-band hmit U <S^ D, Jb.kky/D « 0.62 x 
16V^/U, therefore the critical ti2,c is given by ii2,c ~ T 
and it does not depend on U . This provides an alterna- 
tive interpretation for the [/-dependence of ii2,c in the 
strong-coupling regime found in Ref. 163 . 



C. Isospin-invariance breaking perturbations 

The inter-impurity electron repulsion and the two- 
electron hopping between the impurities represent per- 
turbations that break the isospin SU{2) symmetry of the 
original model, while they preserve both the particle-hole 
symmetry as well as the spin invariance. 



1. Inter-impurity electron repulsion 

The effect of the inter-impurity electron repulsion (in- 
duced by capacitive coupling between the two parallel 
quantum dots) is studied using the Hamiltonian 



-ffdots = y^ Hdot,i + Ui2{ni - l)(n2 - 1), 



(41) 



where it should be noted that (ni — l){n2 — 1) = 4/f/| 
is the longitudinal part of the isospin-isospin exchange 
interaction /i ■12- 

Results in Fig. [TUl show that the inter-impurity repul- 
sion is not an important perturbation as long as U12 < U. 
Finite U12 only modifies the Kondo temperature and the 
temperature T^ of the FO-LM transition, while the be- 
havior of the system remains qualitatively unchanged. 
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Note that Tp is unchanged since U12 equally affects both 
the singlet and the triplet energy. 
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Figure 10: (Color online) Temperature-dependent susceptibil- 
ity of the 2-dot systems for different inter-impurity electron- 
electron repulsion parameters [/12. Circles are the Bethe- 
Ansatz results for the S = 1/2 Kondo model which fit the 
NRG results in the special case f/12 ~ U. 

For U12 > U the electrons can lower their energy by 
forming on-site singlets and the system enters the charge- 
ordering regime^. This behavior bares some resemblance 
to that of the negative- [/ Anderson modetfl which under- 
goes a charge Kondo effect. 

The system behaves in a peculiar way at the transition 
point U12 = U where C/12 and U terms can be combined 
using isospin operators as 



[//2(4(/0'+4(/|)2 



Ui2^im 



2U{F- 



(42) 



We now have an intermediate temperature fixed point 
with a six-fold symmetry of states with /^ = as can be 
deduced from Eq. H42I) and the entropy curve in Fig. ^2 
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Figure 11: (Color online) Temperature-dependent entropy of 
the 2-dot systems for different inter-impurity electron-electron 
repulsion Ui2. 



For two impurities we can define an orbital pseudo-spin 
operator as 



^= 9J2 Yl 4»o■^3dja, 



(43) 



a ij = 1.2 



where <t is the vector of the Pauli matrices. The quan- 
tum dots Hamiltonian -ffdots commutes for U12 — U 
with all three components of the orbital pseudo-spin 
operator; the decoupled impurities thus have orbital 
SU{2)orh symmetry. Furthermore, pseudo-spin O and 
spin S operators commute and the symmetry is larger, 
SU{2)spin <^ SU{2)orh- In fact, the set of three S\ three 
O^ and nine operators S^O^ are the generators of the 
SU{4:) symmetry group of which SU{2)spin (^ SU{2)oih is 
a subgroup. The six degenerate states are the spin triplet, 
orbital singlet and the spin singlet, orbital triplet^^ which 
form a S'[/(4) sextet: 

|5 = i,5, = i,o = o) = |T,T), 
\s = i,s, = o,o = Q) = 1/V2 (I T, i) + I i, T» , 
|5 = 1,5, = -1,0 = 0) = U,i), 

|5 = o,o = i,o, = i) = |Ti,o), 

\s = 0,0 = 1,0^ = 0) = i/V2(| T,i) - I i,T)), 

|5 = 0,O = l,O, = -l) = |0,Ti). 

The states \S = 0,0 = l,Oz = ±1) can be com- 
bined into an isospin triplet jS* = 0,/ = l,/2 = 0) = 
1/V2(| Ti,0) + |0,ti)) and an isospin singlet |S' = 0,/ = 

o) = i/V2(|Ti,o)-|o,n)). 

The coupling of impurities to the leads, however, 
breaks the orbital symmetry. Unlike the model studied in 
Ref. l63l our total Hamiltonian H is not 5'[/(4) symmet- 
ric, so no S'C/(4) Kondo effect is expected. Instead, as the 
temperature decreases the degeneracy first drops from 6 
to 4 and then from 4 to 2 in a S' = 1/2 SU{2) Kondo effect 
(see the fit to the Bcthe-Ansatz result in Fig. llUf) . There 
is a residual two-fold degeneracy in the ground state. To 
understand these results, we applied perturbation the- 
ory (Appendix IXjl which shows that the sextuplet splits 
in the fourth order perturbation in Vk- The spin-triplet 
states and the state 15* = 0, / = 0) form the new four- fold 
degenerate low-energy subset of states, while the states 
\S = 0,1 = l,Iz = 0) and \S = 0,0 = l,Oz = 0) 
have higher energy. The remaining four states can be 
expressed in terms of even and odd molecular-orbitals 
described by operators dl^ = 1/V2 (dl^ + dl^j and 

dl = l/V2(dl-dl). We obtain 



\s = i,s, = i,o = o) = dl^dyo), 

\S=1,S, = 0,0 = 0) = I/V2 (<^4,i + dlAd) |0)' 
\S = 1,S, = -1,0 = 0) = dl^dljO), 

\S = 0,1 = 0) = l/^/2 [di/^^^ - dl^dl^) |0). 

(44) 
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The four remaining states are therefore a product of a 
spin-doublet in the even orbital and a spin-doublet in 
the odd orbital. Due to the symmetry of our problem, 
only the even orbital couples to the leads, while the odd 
orbital is entirely decoupled. The electron in the even 
orbital undergoes S — 1/2 Kondo screening, while the 
unscreened electron in the odd orbital is responsible for 
the residual two-fold degeneracy. 



2. Two-electron hopping 
We consider the Hamiltonian 

2 



-ffdots — 2_^ Hdot,i — T12T, 



(45) 



where T is the two-electron hopping operator that can be 
expressed in terms of the transverse part of the isospin- 
isospin exchange interaction Ii ■ I2: 



f = d\^d\^d2id2ii + dl^dl^diidii 
= I+I2 + I^lt = 2(/f /f + /f /| 



(46) 



This perturbation term is complementary to the one gen- 
erated by Ui2 in Eq. (|^ and studied in the previous 
subsection. Physically, it corresponds to correlated tun- 
neling of electron pairs which can be neglected in the 
applications to problems of transport through parallel 
quantum dots coupled electrostatically as physically real- 
ized in semiconductor heterostructures. Models featuring 
pair-tunneling terms as in Eq. (|46|l may, however, be of 
interest to problems in tunneling through molecules with 
vibrational degrees of freedom, where ground states with 
even number of electrons can be favored due to a pola- 
ronic energy shift^S*Sl. In such cases, the charge trans- 
port is expected to be dominated by the electron-pair 
tunneling^. 

The temperature dependence of the magnetic suscepti- 
bility shown in Fig. ll2l again demonstrates the robustness 
of the S = I state for |Ti2| < U. The behavior of the 
system for negative T12 is similar to the case of the inter- 
impurity repulsion. For T12 = —U we again observe spe- 
cial behavior of the susceptibility curve, characteristic for 
the six-fold degeneracy observed in the previous subsec- 
tion at U12 = U. For positive T12 the system undergoes 
the S = 1 spin Kondo effect up to and including T12 = U. 
The FO-LM transition temperature Tj* and the Kondo 
temperature are largely T12 independent, while the LM- 
FF transition temperature Tp decreases with increasing 

Tl2- 



D. Unequal coupling to the continuum 

We finally study the Hamiltonian that allows for un- 
equal hybridizations F^ = T^po\VJ: p in the following 
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Figure 12: (Color online) Temperature-dependent suscepti- 
bility of the 2-dot system with two-electron hopping between 
the dots T12. 



form: 



with 



H — Hhand + Hdots + 2_^ -^c.i, 



ka- 



(47) 



(48) 



We set 142 = aV,}, i.e. F2 = a^Fi. 

The effective low-temperature Hamiltonian can be now 
written as 

2 

HcS = ffband + S ■ 2_^ JkaSi — ^RKKV^l ' '^2- (49) 
i=l 

with Jk,2 — ce^Jx,! and with the effective RKKY ex- 
change interaction given by a generalisation of Eq. l|^ 



■^RKKY = 0.62UpIJk,iJk,2 = a'jRKKY, (50) 

where Jrkky is the value of RKKY parameter at a = 1 . 
In our attempt to derive the effective Hamiltonian we 
assume that in the temperature regime T < '^rkky ^^^ 
two moments couple into a triplet. Since the two Kondo 
exchange constants Jk.i are now different, we rewrite HcS 
in Eq. (|49|1 in the following form 

^cff = ffband + * ' ( ^ ^ ~ ('^1 + '^2) 

^,fJKS_JK^^g^_g^^\ (51) 



J\ 



cff 
RKKY 



Si ■ S2 



Within the triplet subspace, ^i -I- ^2 is equal to the 
new composite spin 1, which we denote by S, Si — S2 is 
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identically equal to zero, and Si • 5^2 is a constant —1/4. 
As a result, the effective Jk is simply the average of the 
two exchange constants: 



^K.off — 



Jk,1 + Jk,2 



(52) 



Susceptibility curves for different a are shown in 
Fig. ^1 Note that the Kondo temperature determined 
using Eq. jSJ combined with the naive argument given 
in Eq. H52(l fails to describe the actual Kondo scale for 
a < 0.4 as seen from Fig. ^1 This is due to admixture 
of the singlet state, which also renormalizes Jk, even 
though the singlet is separated by Jr^ky ^ "^k from 
the triplet subspace. Note however, that ^rkky i^ ^^^^ 
described by the simple expression given in Eq. (|50|l as 
shown in Fig. 1141 By performing a second-order RG cal- 
culation (see Appendix IB|l . which takes the admixture of 
the singlet state into account, we obtain Tr- as a func- 
tion of a which agrees very well with the NRG results, 
see Fig. [HI 

For extremely small a, ^rkky eventually becomes 
comparable to the Kondo temperature, see Fig. ^^ For 
that reason the ferromagnetic locking-in is destroyed and 
the system behaves as a double 5=1/2 doublet, one of 
which is screened at T^ = Tk{Jk,i) as shown in Fig.lCTp. 



VII. CONCLUSIONS 

We have shown that several magnetic impurities, cou- 
pled to the same Wannier orbital of a conduction electron 
band, experience ferromagnetic RKKY interaction which 
locks local moments in a state of a maximal total spin. 
The multi-impurity Anderson model is at low tempera- 
tures, i.e. for T <Tp, equivalent to a 5 = A^/2 SU{2) 
Kondo model. Using perturbation theory up to the 
fourth order in V we derived an analytical expression for 
>/rkky and tested it against NRG calculations. We have 
also shown that the high-spin state is very robust against 
experimentally relevant perturbations such as particle- 
hole symmetry breaking, on-site energy level splitting, 
inter-impurity capacitive coupling and direct exchange 
interaction. At low temperatures, the ferromagnetically 
locked impurities undergo a collective Kondo cross-over 
in which half of a unit of spin is screened. The Kondo 
temperature in this simple model does not depend on the 
total spin (i.e. on the number of impurities N), while the 
LM-FF temperature Tp is weakly iV-dependent. 

We next list a few most important findings concerning 
the effect of various perturbations to the original two-dot 
system: a) Tp is in the range S < 0.4C/ nearly indepen- 
dent of the deviation from the particle-hole symmetric 
point (5 = 0, b) increasing the difference between on-site 
energies of two dots, 2A, induces a Kosterlitz-Thouless 
type phase transition separating the phase with S — 1/2 
residual spin at low-temperatures from the S = one, 
c) introduction of additional one-electron hopping be- 
tween the impurities induces effective AFM interaction 




0.6 



>0.5 



~S 0.4 



0.3 



0.2, 



1 1 1 1 1 1 1 

— a = 0.1 

— a = 0.05 

— a = 0.01 

— a = 0.005 
. — a = 0.001 

— a = 0.0005 

- — a = 0.0001 J 


1' 


%■- 


a=0.1 // 
S=l Y^ox^Ao //// 


i 


^ a=0.0001 \ 

= 1/2 Kondo \ - 
U/D=0.01,r/U=0.02 


, , 1 ... 1 



25 



-20 



-15 
logio(T/D) 



-10 



Figure 13: (Color online) Temperature-dependent suscepti- 
bility of the 2-dot system with unequal coupling to the leads, 
r2 = Q^Fi. a) The range of q where Tk is decreasing, b) The 
range of a where Tk is increasing again. Circles (squares) are 
BA results for the S = 1 (S = 1/2) Kondo model. The ar- 
rows indicate the evolution of the susceptibility curves as the 
parameter a decreases. 
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Figure 14: (Color online) Comparison of calculated and pre- 
dicted Kondo temperature Tk and effective exchange inter- 
action Jrkky- The calculation of scaling results for Tk is 
described in Appendix IHl 
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Jafm — 'itu/U that does not effect the Kondo tem- 
perature as long as Jafm ^ Jrkky, nevertheless, at 
ti2 — ti2.c it destabilizes the S = \ state. The crit- 
ical value ii2,c ~ T does not depend on [/, d) inter- 
impurity Coulomb interaction U12 leads to a transition 
from the 5=1 Kondo state to the charge ordered state. 
In the 4- fold degenerate intermediate point, reached at 
t/12 = [/, the effective Hamiltonian consists of the ef- 
fective 5=1/2 Kondo model and of a free, decoupled 
5 = 1/2 spin, e) when the two impurities are coupled to 
the leads with different hybridization strengths, second- 
order scaling equations provide a good description of the 
Kondo temperature. 

The properties of our model apply very generally, since 
high-spin states can arise whenever the RKKY interac- 
tion is ferromagnetic, even when the dots are separated 
in spac o*^^'^^ . In addition, it has become possible to study 
Kondo physics in clusters of magnetic atoms on metallic 
surfaccs^-^^. On (111) facets of noble metals such as cop- 
per, bulk electrons coexist with Shockley surface-state 
electronaiSi. Surface-state bands on these surfaces have 
kp ^ 0.1 — 0.2A~^; thus, for nearest and next-nearest 
neighbor adatoms kpR < 1. If hybridization to the sur- 
face band is dominant, small clusters then effectively cou- 
ple to the same Wannier orbital of the surface band and 
the single-channel multi-impurity Anderson model is ap- 
plicable; in the absence of additional inter-impurity inter- 
actions, the spins would then tend to order ferromagneti- 
cally. If hybridization to the bulk band with A:|P ~ 1A~^ 
is also important, the problem must be described using 
a complex two-band multi-channel Hamiltonian. 

Further aspects of the multi-impurity Anderson model 
should be addressed in the future work. Systems of cou- 
pled quantum dots and magnetic impurities on surfaces 
are mainly characterized by measuring their transport 
properties. Conductance can be determined by calcu- 
lating the spectral density functions using the numerical 
renormalization group method. We anticipate that the 
fully screened iV = 1 model will have different temper- 
ature dependence as the under-screened TV > 2 models. 
Since in quantum dots the impurity level S (or e^) can be 
controlled using gate voltages, it should be interesting to 
extend the study to asymmetric multi-impurity models 
for N > 2 where more quantum phase transitions are ex- 
pected in addition to the one already identified for iV = 2 
at (S = Sn- 
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Appendix A: RAYLEIGH-SCHROEDINGER 
PERTURBATION THEORY IN Vk 



Following Ref . |7lJ , we apply Rayleigh-Schrodinger per- 
turbation theory to calculate second and fourth order 
corrections in 14 to the energy of a state \n): 



Ei?^ = E 



'{n\Hc\m){m\Hc\n) 



Ei'^- E ' (Al) 

{n\Hc\rn3){m3\Hc\rn2){m2\Hc\rni){rni\Hc\n) 

(En ~~ E,n3){En " Em2){En — Emi) 



The summation extends over all intermediate states \mi) 
not equal to one of the degenerate ground states. We 
will consider the simplified case of constant 14, i.e. Vk — 

Vk, = V. 



1. RKKY interaction in the two-impurity case 

We study the splitting between the singlet \S) = 
1/V2(| t,i) - U,T)) and the triplet state |r) = 
l/-\/2 (I tii) + l i:?))- The second order corrections are 
^^2) ^ £;^2) = _(5'^ + S2) with 



5,=4y^l E 

k'>k 



L ^f-f U~26 + 2ek' 

k' >kF 



(A2) 



L ^^ U + 25~2ek 

k<.kF 



There is therefore no splitting to this order in V. The 
fourth order corrections are 



4^*) = wf + M^pp + Iy^^ 



E 



(4) 



Wr 



Ph 



(A3) 



T ' 



The authors acknowledge useful discussions with J. 
Mravlje and L. G. Dias da Silva and the financial support 



where the particle-hole (ph), particle-particle (pp) and 
hole-hole (hh) intermediate-state contributions are 
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Wf 



16T/4 1 
U L^ 



E 



£k){5 - ek'){U + ek -Cfe') 



k<~.kF ,k''>ki7 



iU-26 + 2ekY{U + 25- 2ek'Y{ek - £/,') ' 



^ph ^ 32^4^ y 2[/^ + WJEk - efcO + 4 (,5^ + e^ + e^, - ek^k' - Sje^ + ffcQ) 



L2 



-i-i^i. 



k'^kF ,k' 'ykp 

E 

k'^'^kp .k2>kp' 



<-Ti 



i: 



fci <.kF ,k2^kp 



iU-26 + 2ekY{U + 25- 2ek'Y{ek - et') 

2{U - 25){W - 25) + 8{U - ,5)(efc; + e^, ) + Se^^ e^. 
([/-25 + 2efe,)2(C/-25 + 2efe.)2 

2(t/ + 25)iSU + 25) - SjU + 5)iek, + efcj + Sefc.efc, 
{U + 25-2ekJHU + 25-2ek,)^ 

I 



(A4) 



From these expressions we obtain Jrkky = Es — Er- 
in order to evaluate the sums for a flat band with a con- 
stant density of states po = 1/(21?) and the chemical po- 
tential /i = 0, we make formal replacements j^ ^k'>k ~ 



1/2 j^dk' and i^ 



E 



fc<=0 



L 2^fe'>0 ~^ '-/■^Jo ""' '^^^"^ L ^k<kf 

1/2 J .^dk. In Fig. E] we plot the prefactor c in the 
expression for the exchange constant Jrkky = cl6V^ /U 
as a function of U/D. In the wide-band limit, i.e. for 
small U/D, c approaches a constant value of c = 0.616 
irrespective of the value of 5/U. The dependence of c on 
5 for U/D ~ 1 is due to the band-edge effects. 
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Figure 15: The prefactor c in the RKKY exchange constant 
■/rkky — cl6V*/U for a fiat band with po — 1/(2D) as a 
function of U for a range of values of the impurity energy 
level S. 

To determine the temperature Tj at which the RKKY 
interaction becomes fully established, we calculate the 
cut-off dependent Jrkky (.£'), where E is the low-energy 
cut-off for k and k' integrations, i.e. the integrals over 
k and k' become J^dk' and J_^ dk. In Fig. El we plot 

the ratio (,{E) = Jrkky (-E)/J{(kky: where Jrkky = 
Jrkky (£' ~^ 0). The ratio £,{E) reaches an (arbitrarily 
chosen) value of 0.9 at E/U ~ 0.02. This value of E 
roughly defines Tj below which the RKKY is fully de- 



veloped. For small enough V (i.e. F), the value of Tj 
is positioned between Tj* (free-orbital to local-moment 
transition temperature) and Tp, the temperature of fer- 
romagnetic ordering of spins, given by Tp — Jrkky//^' 
where /? is a constant of the order one. For larger V , 
however, Jrkky (J") does not reach its limiting value at 
the temperature where the spins start to order. In this 
case we obtain the ordering temperature Tp numerically 
as the solution of the implicit equation 



Tp — Jrkky (Jf)//3- 



(A5) 



An approximate fit to S,{E) in the wide-band limit is 
^{E) = 1/(1 + xE/U) with X = 12.2. We then obtain 



a solution for Tp in closed form: 



^, _ v/l + 4x//3(J|JKKY/t^) 
^ 2x/U 



- 1 



2x/U 

, -^RKKY I 1 „ ^ '^RKKY i /O ( '^RKKY ] ] 

' P \ p u ^ \ u ) )■ 



(A6) 
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Figure 16: Ratio S,{E) = Jrkky (i?)/ Jrkky of the running 
RKKY coupling constant at energy E over its value in the 
E —* Q limit. The dashed line is an approximate fit to a 
simple rational function (,{E) — 1/(1 + xE/U). 
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2. Six-fold symmetric U12 = U case 

We study the splitting between the singlet, the 
triplet (same as above) and the "exciton" states 
|J = 0) = l/A/2(|n,0)-|0,n)) and |/ = 1) = 
l/\/2 (I t j, 0) + |0, ti))- Second order corrections are all 



equal: 



p(2) _ p(2) _ 



E 



(2) 
1=0 



= E 



(2) 
7=1 



-(^1 + ^2) 



where 5*1 and S2 are the same as in the previously treated 
C/12 = case. There is again no splitting to second order 



in V . The fourth order corrections are 



where 



Ef 


= Wf + T4^|P 


p(4) 

H/rp 


= M^P^ 


p(4) 
^7=0 


p(4) 


p(4) 


= 4^) 



w, 



hh 



(A7) 



J 



W. 
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k'Ckp ,k' '^kf 
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W^"" = 32y4i_ ^ 
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W^ + 6£/(efc - efcO + 4 (^^ + e^ + e^, - ekeu' - 5{^k + £fc')) 
{U~25 + 2ekf{U + 25-2ekf(ek-ek') 

2(C/-25 + efe;+efc^)2 



WPP = lfiF4j_ Y^ 

' ^\'~..v^, {U-25 + 2ek0^{U-25 + 2ek0H2U-25 + ek'^+ek0' 

ki >kF ,k2 >kF \. z 
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hh 



IQV^ 
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2{U + 25-€k,+ek,f 
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^. ^.^^i (C/ + 25-2efeJ2([/ + 2(5-2efeJ(2t/ + 2(5-efe, +efej- 
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The triplet is degenerate with the / = state, while 
the singlet and the 1 = 1 state are higher in energy, as 
determined by performing the integrations (results not 
shown) . 



Appendix B: SCALING EQUATIONS TO SECOND 
ORDER IN J 



We consider an effective Hamiltonian of the form 



ka 



/ , •Jmm'-^inm'C^^Ck'a' , 



(Bl) 



vam'kk'1717' 



where Xmm' — \'m){m'\ are the Hubbard operators and 
^rmn' ^'^^ generalized exchange constants. 
We writeS 



Hn + H12 {E — H22) H21 



Hio {E — Hqo) Efoi 



Tpi = Eipi, 



(B2) 



where subspaces 2 corresponds to states with one elec- 
tron in the upper \5D\ edge of the conduction band, 
corresponds to states with one hole in the lower \5D\ 
edge of the band, and 1 corresponds to states with no 
excitations in the edges that are being traced-over. Fur- 
thermore, Hij = PiHPj, where Pi are projectors to the 
corresponding subspaces i. 



To second order, the coupling constant are changed by 



SJZa' = Po\SD\Y, 



E-D + Ek-En-Ho 



jTCr jT(T 



-pomj: ^_^_J_ 



E-D-Ek'- E„ - H, 



jCr T J(7T 



(B3) 



We apply these results to the effective low-temperature 
Kondo Hamiltonian 

-ffcff = ifband + JlS-5l + J2S-52- JgW (-^l ' -^2 - 1/4) . 

(B4) 
Introducing spin-1 operator S defined by the following 
Hubbard operator expressions: Sz = X^-^ — X^^ S^ = 
%/2 (X|o + ^oi) and S^ = {S+)\ we obtain 

H = i?band + Js ■ S + Jy{KKyXsS 

+ A{sz{Xos + Xso) + s+{X^s-Xs^) (B5) 

+ s-{Xsi-X^s)), 

where index S denotes the singlet state and we have 

J1 + J2 



J 

A 



2 
Ji- J2 



Jo(l + a')/2, 
Jo(l-a')/2. 



(B6) 
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Equations (|B3p reduce to two equations for J and A 
"j2 A2 



5J = pa\5D 
SA 



u u ^ '^RKKY 



(B7) 



2«MaX 



from which ensue the foUowing scahng equations 



where I = logD. The initial bandwidth D is the effective 
bandwidth I?off = 0.182C/ for the Anderson model and we 
take J {I = logZ?cff) = J and A(Z = log Deft) = A with J 
and A taken from Eq. IIB6|l . We integrate the equations 
numerically until J starts to diverge. The corresponding 
cut-off D defines the Kondo temperature. The results 
are shown in Fig. 1141 The scaling approach reproduces 
our NRG results very well. 



dJ _ j2 ^^D 



dA ^ , ~ 



RKKY 
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